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Abstract 

Firstly, we reduce the long-standing problem of ascertaining the Hilbert-Schmidt probability that 
a generic pair of qubits is separable to that of determining the specific nature of a one-dimensional 
(separability) function of the radial coordinate (r) of the unit ball in 15-dimensional Euclidean 
space, and similarly for a generic pair of rebits, using the 9-dimensional unit ball. Separability 
probabilities, could, then, be directly obtained by integrating the products of these functions 
(which we numerically estimate and plot) with jacobian factors of r m over r £ [0, 1], with m = 17 
for the two-rebit case, and m = 29 in the two-qubit instance. Secondly, we repeat the analyses, 
but for the replacement of r as the free variable, by the azimuthal angle 4> E [0, 27r]-with the 
associated jacobian factors now being, trivially, unity. So, the separability probabilities, then, 
become simply the areas under the curves generated. For our analyses, we employ an interesting 
Cholesky-decomposition parameterization of the 4x4 density matrices. Thirdly, in an exploratory 
investigation, we examine the Hilbert-Schmidt separability probability question using the three- 
dimensional visualization (cube/tetrahedron/octahedron) of two qubits associated with Avron, 
Bisker and Kenneth, and Leinaas, Myrheim and Ovrum, among others. We find-in its two-rebit 
counterpart-that the Hilbert-Schmidt separability probability-the ratio of the measure assigned 
to the octahedron to that for the tetrahedron-is vanishingly small, in apparent conformity with 
observations of Caves, Fuchs and Rungta regarding real quantum mechanics. However, the ratio of 
the measure assigned to entangled states to that for potential entanglement witnesses (represented 
by points in the cube) is small, ~ 0.005, but seemingly finite. We, then, begin an investigation of 
the full 15-dimensional two-qubit form of the question. 
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I. INTRODUCTION 



A 2006 monograph entitled "Geometry of Quantum States: An Introduction to Quan- 
tum Entanglement" pQ [2] had received by September 2011 (according to scholar. google) 
more than 460 citations in the scientific literature. We aspire, in this communication, to 
contribute-in two distinct sets of analyses (sees, [n] and III )-to this general area of strong 
current interest. As in a number of previous studies [3H9], we will be interested in deter- 
mining (univariate) "separability functions", and with their use, "separability probabilities" 
[TU] . However, the analytical (Cholesky-decomposition) parameterization of density matri- 
ces which we, firstly, employ in secjllj is different than the (Euler-angle, Bloore/correlation) 
parameterizations earlier utilized, and our results here reflect different aspects of the one and 
the same, fundamental underlying, still not fully resolved, separability probability problem. 



II. CHOLESKY-DECOMPOSITION PARAMETERIZATION ANALYSIS 

We begin by noting that in three preprints in the past decade, Ramakrishna and various 
colleagues have examined a number of parameterizations of nonnegative-definite matrices 
(which include, of course, among them, the density matrices of quantum theory-bearing 
the additional requirement beyond nonnegative-definiteness of unit trace) [TTHTB"] . One of 
these parameterizations-that will be the focus of our study here-is based on the well-known 
Cholesky-decomposition (cf. [TlTflB"] ). which is of important use in the numerical solution 
of linear equations, among other areas. 

We can represent a density matrix p, using the decomposition, as the product of a lower- 
triangular matrix T = HT^H-so that T^j = 0, j > fc-with nonnegative entries (r^ > 0) on 
the diagonal, and the transpose of its complex conjugate. Thus, 

p = rrt. (i) 

It can be seen that the determinant of an n x n density matrix p, having imposed the unit 
trace requirement p nn = 1 — ^ZiPkk, takes the form (cf. (T2J p. 8] for a different form of 
determinantal expression; also [TFJ sees. D.2 and D.3]), 
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Further, 

Pnn = (1 — ^k=l^kk ~ ^k=l^ k j<k\^ kj\ 2 ) (3) 

Thus, if the second factor in ^ is nonnegative, p nn -which is at least as large as this factor- 
must also be nonnegative, and all the requirements that p be a density matrix are fully 
satisfied. One readily perceives, then, that for the case n = 4, the two-qubit states (p) can 
be put into one-to-one correspondence with the points-imposing the triad of restrictions 
Tfcfc > (k = l,2,3)-of one-eighth of the unit ball (centered at the origin, of course) in 
15- dimensional Euclidean space. Further, for the two-rebit states-that is those assigned 
4x4 density matrices, the entries of which are restricted to real values [T5H2U] — -the one- 
to-one correspondence is with the points of an analogous one-eighth part (a "hemidemisemi- 
hypersphere" ) of the unit ball in 9-dimensional space. 

The jacobian for the transformation to the lower-triangular T parameters is 

j dCtwo-rebit = ^-^^22^33 (4) 

in the two-rebit case, and 

JClCtwo-qubit = ^ 22^ 33 (5) 

in the two-qubit case. 

So we see, it is quite natural-given the Cholesky-decomposition parameterization and 
the sum-of- squares nature of the second factor in (|2]) to convert the T's to hyperspherical 
coordinates (r, Q\, 62, ■ ■ •) (cf. [15, eq. (A. 2)]). For the radial coordinate (r) we can then 
write (cf. @), 

r = ^=ln k + ^=^U\ T ^ 2 ( 6 ) 
in these two differing (n = 9,15) dimensional settings. (Of course, |rfcj| 2 = T\- in the 
two-rebit case.) 

We should, importantly, point out that the origin of the unit balls in both the two-qubit 
and two-rebit cases does not correspond, as one might initially expect, to the classical mixed 
state-having a diagonal density matrix with all entries equal to |-but to a [pure] diagonal 
density matrix with one entry equal to 1 and the other three to 0. By specifically choosing 
to set p 44 = 1 — /On — P22 ~ P33, we have in our analyses that this unit diagonal-entry lies in 
the (4,4)-position. 
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Our principal objectives in this section of the study are to determine the relative 
prevalence-with respect to the Hilbert-Schmidt measure [I, sec. 14.3] [211 [22]~of separa- 
ble vs. nonseparable/entangled states as: (1) the radial coordinate r of the unit ball varies 



between (corresponding to the origin of the unit ball) and 1 (its boundary) (sec. II A); 



and (2) the azimuthal coordinate varies between and 20 (sec. II B). From the results, 
we will be able to directly (through univariate numerical integrations) obtain estimates of 
the probabilities of separability [10]. We also aspire to discerning if the associated curves 
(separability functions) appear to assume any special functional form-in which case, the 
separability probabilities themselves might be perceptible. 



A. Radial Analyses 

We can now make an important observation. The requirement that the partial transpose 
of p have a nonnegative determinant-and, thus, corresponds to a separable state |23j-takes 
the general (quadratic) form (in both the two-qubit and two-rebit cases), 

f(e l ,9 2 ,...)-g(e 1 ,e 2 ,..y>o. (7) 

Here, of course, the (apparently quite complicated) functions / and g of the angular (non- 
radial) coordinates are specific to the two-qubit and two-rebit scenarios. (Obviously, for 
/ = g — 1, (J7| simply reduces back to the hyperspheres [unit balls] in question-one-eighth 
parts of which are, as already noted, the loci of the associated feasible 4x4 density matrices.) 

If we enforce the separability requirement ([7]), while numerically integrating the corre- 
sponding jacobians over the angular coordinates (fourteen in number in the two-qubit case, 
and eight in the two-rebit case), and normalize the results by the corresponding Hilbert- 
Schmidt volumes [2T]-that is, 96 yg 80 (two-rebit volume) and 10897 2864ooo (two-qubit volume)- 
we obtain the curves (radial profiles) displayed in Figs. [l]and[2j 

The integrands used in this process are the products of the indicated jacobians (Q,([5j) 
with the jacobians for the transformations from the lower-triangular parameters T^j to hy- 
perspherical coordinates. (The exponent of r in the resultant composite two-rebit jacobian 
is 17=9+8, while it is 29=15+14 in the two-qubit counterpart. The integers 8 and 14 are 
contributed by the standard jacobians for the transformations from cartesian to hyperspher- 
ical coordinates , and the integers 9 = 4 + 3 + 2 and 15 = 7 + 5 + 3 from the j acobians ( Q , (pj ) 
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for the transformation to Cholesky-decomposition variables.) We employ numerical [Monte 
Carlo] integration-with 10,600,000 fourteen-dimensional points in the two-qubit case, and 
1,100,000 eight-dimensional points in the two-rebit case. (More points were used in the for- 
mer case, since due to the "curse-of-dimensionality", convergence to a smooth separability 
function was decidedly slower there. The use of quasi-Monte Carlo procedures (2U 125]- 
which we will use for Figs. [7] and [9]-might also have been somewhat more computationally 
effective.) 

For each individual randomly-generated point, we evaluate whether the separability con- 
straint <Q is satisfied or not at each of 10,001 equally-spaced values of r 6 [0, 1]. (We have 
also attempted to determine-by enforcing, through symbolic means, the separability con- 
straint ([7]) with r = 0-the exact theoretical values of the y-intercepts in the two figures, but 
have so far been unsuccessful in this quest. Clearly, the enforcement in symbolic terms of ^ 
with r free to vary-that is, the separability probability question in its full generality-seems 
even more computationally-challenging.) 

The two-rebit curve in Fig. [TJ if scaled to equal 1 at r = 0, can be very well-fitted, we 
found, by a curve (Figs. M and p} 



1 f r t 2 (l-t)- 3 / 4 1 

I r (r, 3, -) = 1 - \^ - = h- r (l - r, -, 3). (8) 

V 4 ; /V(l-t)- 3 / 4 4' ; y J 



Here I r (r, a, b) is a regularized incomplete beta function 

Let us jointly plot (Fig. [5]) the two-qubit curve (Fig. [2]), along with the |-power of the 
two-rebit curve (Fig. [T|, after scaling both curves to assume the value 1 at r = 0. The 
rather close fit here is somewhat puzzling, as earlier studies of ours (reviewed, in detail, 
immediately below) of other "separability functions" would seem to suggest-in line with the 
"Dyson-index ansatz" of random matrix theory-that a power of 2 would be the appropriate 
one HE]. 



1. Relations to earlier analyses 

The curves in our various figures might be termed "separability functions" of the radial 
coordinate r. From such a perspective, they fulfill a role strongly analogous to the (also 
univariate) separability functions estimated using density-matrix parameterizations (Eu- 
ler angle, Bloore/correlation)-other than the Cholesky-decomposition one used above-in a 
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FIG. 1: The relative prevalence-with respect to the Hilbert-Schmidt measure-of two-rebit separa- 
bility, as a function of the radial coordinate r of the unit ball in 8-dimensional Euclidean space. 
The two-rebit separability probability estimate-0.465885-is obtained by integrating the product of 
this curve and r 17 over r G [0, 1]. The exact probability conjectured in [7J was ^ ~ 0.470588. 




i.n 



FIG. 2: The relative prevalence-with respect to the Hilbert-Schmidt measure-of two-qubit separa- 
bility, as a function of the radial coordinate r of the unit ball in 15-dimensional Euclidean space. 
The two-qubit separability probability estimate-0.209417-is obtained by integrating the product 
of this curve and r 29 over r £ [0, 1]. The exact probability conjectured in [7J was ^ ~ 0.242424. 



number of related studies of ours. [3H9], 

When employing the Bloore/correlation framework, the variable of interest (the counter- 
part to the radial coordinate r above) was the diagonal-entry ratio (or its square root or 
logarithm) of v = • When relying upon Euler angles [27] , the variable of interest was the 
maximal concurrence 



that is, C max = max (Ai — A 3 — 2a/A 2 A 4 ), where the A's are 
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FIG. 3: The (two-rebit) curve in Fig. [TJ scaled to equal 1 at r = 0, plotted along with the closely- 
fitting function (fsl) 




-0.005 - 



-0.010 



FIG. 4: The (two-rebit) curve in Fig. [TJ scaled to equal 1 at r = 0, minus the fitted function ^ 




FIG. 5: A joint plot of the two-qubit (rougher) curve (Fig. [2]) and the |-power of the two-rebit 
curve (Fig. [TJ-after scaling them both to equal 1 at r = 0. 
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the eigenvalues of p, ordered so that Ai > A2 > A3 > A4. In the former (Bloore/correlation) 
set of analyses, the two-qubit separability function gave evidence of being proportional to 
the square of the two-rebit separability function. In the latter (Euler-angle) set of analy- 
ses, there were-quite remarkably and still unexplainedly-pronounced jump discontinuities-in 
both the two-rebit and two-qubit cases-in the associated separability functions at C max = | . 
For C m ax £ [§, 1], but not for C max G [0, |], the two functions appeared to adhere-as in the 
Bloore CcLS6, clS just noted-to the Dyson- index ansatz (/3 = 1,2,4) of random matrix theory. 

In the interest of completeness, in particular in light of the (|)-power relationship appar- 
ently adhered to in Fig. [51 let us review the results given in jT] pertaining to separability 



functions parameterized (in the Bloore/correlation framework) by v = mP44 . The two- 

^ y 1 P22P33 

rebit jacobian based on the transformation to the v variable of the Hilbert-Schmidt volume 
element (P11P22P33P44) 3 / 2 [22] took the form, 

z/ 3/2 (12 {y{y + 2) {y 2 + 14^ + 8) + 1) log (yg - 5 (5u 4 + 32z/ 3 - 32u - 5)) 

JrealKV) ~ 3780(^-1)9 ' ^ 

while the separability function-based on numerical evidence-was surmised to be 

l/-,2)= 1 -(3-u) v ^. (10) 

One can perform the integration 

f 1 111 

2 / J rea i(v)lJv,-,2)dv = = — — - — , (11) 

J y 1 y '2' ; 151200 2 5 -3 3 -5 2 -7' y ' 

and then obtain the conjecture made in [7] for the Hilbert-Schmidt two-rebit separability 
probability, that is 

8 2419200 f 1 „ . w . 1 „, , 

= J q JrealWAv, ^2)dv . (12) 

Also, the further conjecture was made for the Hilbert-Schmidt two-qubit separability prob- 
ability (note-pursuant to the Dyson-index ansatz-involving the square of the function), 

8 48432384000 f 1 „ , _ , 1 



33 17 



f 1 

/ Jcomplex{.v)I u (v,-,2) 2 dv. (13) 
JO 1 



Here J CO mpiex(y) is the complex counterpart, the two-qubit jacobian based on the transfor- 
mation to the v variable of the Hilbert-Schmidt volume element (P11P22P33P44) 3 |22j-with a 
somewhat more intricate expression- to (g. 
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In [20], it was found in the generic two-rebit case that the determinants \p\ and \p 



PT\ 



("PT" denoting partial transposition)-see (16), (17) below-were orthogonal over the associ- 
ated nine-dimensional space with respect to Hilbert- Schmidt measure. Now, having recog- 
nized (sec. [Tl]) that this space can-using the Cholesky-decomposition parameterization-be 
viewed as a one-eighth part of the nine-dimensional unit ball (cf. eq. (A,2)]), it may be 
more readily possible to classify these multivariate orthogonal polynomials (MOPS) [30l I3T] 
into standard forms (cf. [32]). 



B. Azimuthal Analyses 

Now, let us repeat our two-qubit and two-rebit analyses above, but now using-instead 
of the radial coordinate r-as our free variable, the azimuthal angle (0) in each of the two 
hyperspherical coordinate systems. (The counterparts of the separability requirement (J7]) 
now take the somewhat more complicated forms of quadratic polynomials in sin and cos 0.) 
Then, is the only angle with a range [0, 2ir] (rather than [0, 7r]), and also the only angle 
absent from the associated jacobians over which the multidimensional integrations must be 
performed. 

Therefore, to obtain separability probability estimates from the associated "azimuthal 
profiles" we will only need to integrate the profiles of = ^- over G [0, 1] without 
multiplication by any (non-trivial) jacobian factor at all (which were r 17 and r 29 in the 
radial coordinate case). (For each sample point generated, we now evaluate the determinant 
at each of 1,001 equally-spaced values of 0, rather than the 10,001 equally-spaced values of r 
used in the earlier set of analyses. For the two-rebit case, 22,825,000 eight-dimensonal Monte 
Carlo-generated points were employed, while for the two-qubit case, 17,460,000 fourteen- 
dimensional quasi-Monte Carlo points were used-given that simple Monte Carlo results in 
this latter case seemed highly unstable in nature. Nevertheless, even with this change in 
computational strategy, considerable instability appears to remain.) 

Thus, the separability probability estimates would simply be the areas under the 0- 
parameterized separability functions (Figs. [6] and [7])-if the horizontal axes were drawn to 
intercept at F((j>) = 0. (At such a resolution, the two curves look largely flat in nature.) In 
the two-rebit case, there appears to be oscillatory behavior symmetrical (perhaps sinusoidal 
in nature) around = | (or = it) (Fig. 8). In fact, it seems possible to consider Fig. 6 as 
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FIG. 6: The relative prevalence-with respect to the Hilbert-Schmidt measure-of two-rebit separa- 
bility, as a function of the azimuthal coordinate </> = of the unit ball in 8-dimensional Euclidean 
space. The two-rebit separability probability estimate-0.433082-is obtained by integrating this 
curve (the jacobian factor is trivially unity) over G [0,1]. The exact probability conjectured in 
was A ps 0.470588. 
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FIG. 7: The relative prevalence-with respect to the Hilbert-Schmidt measure-of two-qubit sepa- 
rability, as a function of the azimuthal coordinate (ft of the unit ball in 15-dimensional Euclidean 
space. The two-qubit separability probability estimate-0.289656-is obtained by integrating this 
curve (the jacobian factor is trivially unity) over (ft G [0,1]. The exact probability conjectured in 
was | « 0.242424. 

consisting of a single flat baseline, with the oscillations around the baseline simply averaging 
out to zero. In such a case, it would be irrelevant to the chief goal of determining the HS 
separability probability, what the specific amplitude of these (zero-average) oscillations about 
the determinative baseline is. Fig. [7] serves as the two-qubit analogue of Fig. [6j 
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FIG. 8: A fit to the curve (Fig. ^ by the oscillatory function 0.433082 + 0.0244966 cos 4?r0 

We have, above, used as free variables, the radial and azimuthal coordinates. Of course, 
we might, similarly examine any or all of the remaining (non-azimuthal, polar) angles (13 
[two-qubit] or 7 [two-rebit] ones). In the four distinct (two radial in nature, and two, 
azimuthal, as well as two two-rebit and two two-qubit) Monte-Carlo and quasi-Monte Carlo 
analyses, reported above, we would certainly have wished to have been able to employ even 
more points than we have, in the face of computational limitations, been able to. (In any 
case, the computational resources employed-a number of months on MacMini machines-have 
been quite extensive.) 

III. THREE-DIMENSIONAL VISUALIZATION-BASED ANALYSES 

Another interesting, related topic we have investigated is the possibility of ad- 
dressing the Hilbert- Schmidt separability probability problem within the "visualization" 
(cube/tetrahedron/octahedron) framework of Avron, Bisker and Kenneth [33j EI] (cf. [35l 
eq. (1)] [3ni Fig. 4], also [3TH31?] ). ("The octahedron represents the equivalence class of 
separable states. The set of points that lie outside the octahedron but inside the tetrahe- 
dron represent the equivalence class of entangled states. The set of points that lie outside 
the tetrahedron but inside the cube represent entanglement witnesses. The vertices of the 
tetrahedron represent the equivalence class of pure states. Points related by the tetrahedral 
symmetry represent the same equivalence class" [34~].) 

In this regard, we present Fig. |9j based on 128,000 six-dimensional quasi-Monte Carlo 
points (for each such point, we compute the Hilbert-Schmidt measure on a 51 x 51 uniform 
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grid on the d\ — (^-plane-omitting the G?2-axis, since the measure proves to be independent 
of d 2 ) parameterizing the tensor product A <g> B of equation (1) of [35] 

,_ {A®B)p{A®By 
9 Tr[(A®B)p(A®Byy { ] 

where {p\. is a Pauli matrix), and 

P=\(l+ Zl =1 d i( r k <g> <7 fc ) , dj G [-1, 1]. (15) 

(The 2x2 three-parameter [unit determinant] matrices A and B belong to the special linear 
group SL(2,R). The group itself has infinite volume [101 App. 2].) 

The octahedron, tetrahedron and cube are, of course, three-dimensional objects. But 
we have found that the Hilbert-Schmidt measure-in the generic two-rebit case-is constant 
across one of the three-dimensions (d 2 in the notation of sec. IV,B]). (Certainly, it 



is rather obvious-in partial retrospect-that if d 2 ^ in (15), then p must have complex 
entries, and thus does not represent a two-rebit system.) So, we conveniently-for visual 
purposes-require only two dimensions (di,d 3 ) for the plot. For our parameterization of 
elements of the special linear group SL(2, R) we employ a certain formula jHJ eqs. (1-37) 
and (1.38), p. 85]-drawn from Bargmann |42j-because for numerical purposes, all three of 
the parameters have finite range. 

The exceptional feature here-in contrast to our other investigations of "separability 
functions" -is that there is only one germane function-not two distinct ones (one for all 
states, and one just for separable states)-over the smaller number of parameters that is the 
domain of the separability function. 



In this regard, we note that (cf. (14)) 



■ / 1 _ (di — d 2 — d 3 — 1) (c?i + d 2 - d 3 + l) (d 1 - d 2 + d 3 + 1) (d x + d 2 + d 3 - 1) 

|P 1 ~ 256(Tr[(A®5)p(A®5)t])4 t b J 

and, for the partial transpose of p , the slightly different expression (note the sign of 1 in 
the last factor of the two determinantal expressions) 

PT , _ (di —d 2 — d 3 + l) (di + d 2 — G? 3 — 1) (di - d 2 + d 3 - 1) {d\ + d 2 + d 3 + 1)) 



256(Tr[(A® B)p(A® 5)t]) 4 

(17) 



(As already noted, these two determinants have been found to be mutually orthogonal with 
respect to Hilbert-Schmidt measure [20].) Since the denominators of these last two deter- 
minantal expressions are necessarily positive, we see that the positivity of the determinants 
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FIG. 9: The estimated Hilbert-Schmidt measure over the d\ — c^-plane (cf. Fig. 4 of [36] and 
Figs. 1 of both |33} [34J) of the three-dimensional cube of equivalence classes of two rebits. The 
measure is constant across the omitted coordinate d^- The associated separability probability-in 
conformity to real quantum mechanics [18]— is zero. 

themselves depends only upon the numerators-which, importantly and exceptionally, are 
functions only of the three cfj's, and not of any of the six SL(2, R) parameters (for the 2x2 
matrices A and B). 

A. Significance of real quantum mechanics 

Our estimate-based on the results used in Fig. [9]-of the ratio of the Hilbert-Schmidt 
measure assigned to the entire cube representing the equivalence classes of potential entan- 
glement witnessess to the tetrahedron of all possible (entangled and non-entangled) states 
is 0.00584127. However, the estimates of the separability probability (the ratio of the Rie- 
mannian volume of the octahedron to that of the tetrahedron)-to our initial considerable 
surprise-was not on the order of 0.45 as we have repeatedly observed, as noted above, in 
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prior studies. We now found this probability to be vanishingly small. Nevertheless, these 
two disparate outcomes do seem reconcilable using a perspective provided by Caves, Fuchs 
and Rungta [THJ (cf. 031). 

The crucial observation of these three authors-in this regard-is given in the penultimate 
paragraph of [18] : 



Another interesting fact is how the regions of entangled vs. separable states 
within the full set of quantum states differ in going from real to complex quantum 
mechanics. In the complex theory, the maximally mixed state p = \l± of two 
qubits is surrounded by an open set of separable states .... In the real theory, 
however, the states 



AB 



-(i <S> I + a(a y ® aj,)) , a e [0, 1] 



demonstrate that there are entangled states arbitrarily close to the maximally 
mixed state. This conclusion also follows from the fact that the condition for 
real (emphasis added) separability, tr(p AB a y ® a y ) = 0, implies that the space 
of real (emphasis added) separable states is an 8-dimensional submanifold of the 
9-dimensional manifold of all states. 

Thus, it can be seen that our earlier estimates of a Hilbert-Schmidt separability proba- 
bility of a pair of rebits on the order of 0.45 have been (quite legitimately, it would seem) 
obtained in the framework of complex quantum mechanics, while the vanishingly small prob- 
ability observed in the analysis here, relying upon the two-qubit visualization scheme p3H36] , 
is a result of real quantum mechanics-since an 8-dimensional submanifold must necessarily 
have zero 9-dimensional volume. 

Of course, it would be of obvious interest-but computationally quite challenging-to ad- 
dress these measure-theoretic questions in the full two-qubit (rather than the somewhat 
degenerate two-rebit specialization) framework. We are presently investigating this matter. 
In fact, let us present here some very preliminary results in this regard. We have employed 
114,000 15-dimensional quasi-Monte Carlo points, parameterizing the two (six-dimensional) 



general linear groups employed in (14), as well as the three variables di,d2,d%. (We allow the 
entries of the 2x2 matrices A and B to vary uniformly between -500 and 500.) We, further, 
bin the values of dd^-, d% G [—1, 1] generated into 343 = 7 3 uniformaly-sized bins. For each of 
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FIG. 10: Estimated Hilbert-Schmidt measure over the d,2 — c^-plane for two-qubit systems 

these 15-dimensional points, we calculate the corresponding 15 x 15 jacobian determinant. 
(An interesting, importrant question that arises is under what conditions this determinant 
is negative (cf. |44j).) The estimate of the Hilbert-Schmidt two-qubit separability probabil- 
ity obtained is 0.07954. Summing the bins over each of the three coordinates c?j, in turn, 
and interpolating the results over the [—1, l] 2 square, we obtain the three two-dimensional 



(marginal) surfaces depicted in Figs. 10, 11 and 12 We will continue to add quasi-Monte 
Carlo points to these analyses, as well as conduct parallel analyses with more refined bin 
sizes than above. 



Acknowledgments 

I would like to express appreciation to the Kavli Institute for Theoretical Physics (KITP) 
for computational support in this research, as well as to Joseph Avron and Oded Kenneth, 
and Sergio Cacciatori for relevant communications. 



[1] I. Bengtsson and K. Zyczkowski, Geometry of Quantum States (Cambridge, Cambridge, 2006). 

16 



FIG. 11: Estimated Hilbert-Schmidt measure over the d\ — c^-plane for two-qubit systems 

[2] P. B. Slater, MathSciNet p. MR2230995 (2007k:81001) (2007). 

[3] P. B. Slater, J. Phys. A 43, 195302 (2010). 

[4] P. B. Slater, J. Phys. A 42, 465305 (2009). 

[5] P. B. Slater, J. Phys. A 41, 505303 (2008). 

[6] P. B. Slater, J. Geom. Phys. 59, 17 (2009). 

[7] P. B. Slater, J. Phys. A 40, 14279 (2007). 

[8] P. B. Slater, J. Geom. Phys. 58, 1101 (2008). 

[9] P. B. Slater, Phys. Rev. A 75, 032326 (2007). 
[10] K. Zyczkowski, P. Horodecki, A. Sanpera, and M. Lewenstein, Phys. Rev. A 58, 883 (1998). 
[11] T. Constantinescu and V. Ramakrishna, quant-ph/0306167. 
[12] M. C. Tseng, H. Zhou, and V. Ramakrishna, quant-ph/0610020. 
[13] M. C. Tseng and V. Ramakrishna, quant-ph/0610021. 

[14] K. Banaszek, G. M. D'Ariano, M. G. A. Paris, and M. F. Sacchi, Phys. Rev. A 61, 010304 
(1999). 

[15] J. Daboul, Nucl. Phys. B4, 180 (1967). 

[16] F. Verstraete, K. Audenaert, J. Dehaene, and B. D. Moor, J. Phys. A 34, 10327 (2001). 



17 



FIG. 12: Estimated Hilbert-Schmidt measure over the d\ — c^-plane for two-qubit systems 
[17] P. B. Slater, arXiv: 1109.2560. 

[18] C. M. Caves, C. A. Fuchs, and P. Rungta, Found. Phys. Letts. 14, 199 (2001). 
[19] J. Batle, A. R. Plastino, M. Casas, and A. Plastino, Opt. Spect. 94, 1562 (2003). 
[20] P. B. Slater, arXiv: 1007.4805. 

[21] K. Zyczkowski and H.-J. Sommers, J. Phys. A 36, 10115 (2003). 
[22] A. Andai, J. Phys. A 39, 13641 (2006). 

[23] R. Augusiak, R. Horodecki, and M. Demianowicz, Phys. Rev. 77, 030301 (R) (2008). 
[24] G. Okten, MATHEMATICA in Educ. Res. 8, 52 (1999). 

[25] H. Faure and S. Tezuka, in Monte Carlo and Quasi-Monte Carlo Methods 2000 (Hong Kong), 
edited by K. T. Tang, F. J. Hickernell, and H. Niederreiter (Springer, Berlin, 2002), p. 242. 

[26] A. K. Gupta and S. Nadarajah, Handbook of Beta Distribution and Its Applications (Marcel 
Dekker, New York, 2004). 

[27] T. Tilma, M. Byrd, and E. C. G. Sudarshan, J. Phys. A 35, 10445 (2002). 

[28] F. Verstraete, K. Audenaert, and B. DeMoor, Phys. Rev. A 64, 012316 (2001). 

[29] R. Hildebrand, J. Math. Phys. 48, 102108 (2007). 

[30] I. Dumitriu, A. Edelman, and G. Shuman, J. Symb. Comp. 42, 587 (2007). 



18 



[31] C. Dunkl and Y. Xu, Orthogonal Polynomials of Several Variables (Cambridge, Cambridge, 
2001). 

[32] R. C. Griffiths and D. Spano, Bernoulli 17, 1095 (2011). 

[33] J. E. Avron, G. Bisker, and O. Kenneth, J. Math. Phys. 48, 102107 (2007). 

[34] J. E. Avron and O. Kenneth, Ann. Phys. 324, 470 (2009). 

[35] F. Verstraete, J. Dehaene, and B. D. Moor, Phys. Rev. A 64, 010101(R) (2001). 

[36] J. M. Leinaas, J. Myrheim, and E. Ovrum, Phys. Rev. A. 74, 012313 (2006). 

[37] S. Szarek, E. Werner, and K. Zyczkowski, J. Math. Phys. 49, 032113 (2008). 

[38] S. Szarek, E. Werner, and K. Zyczkowski, J. Phys. A 44, 045303 (2011). 

[39] M. D. Lang and C. M. Caves, arXiv:1006.2775. 

[40] B. Simon and M. Taylor, Commun. Math. Phys. 101, 1 (1985). 

[41] W. Miller, Jr., Symmetry and Separation of Variables (Addison- Wesley, Reading, 1977). 
[42] V. Bargmann, Ann. of Math. 48, 568 (1947). 

[43] J. Batle, A. R. Plastino, M. Casas, and A. Plastino, Phys. Lett. A 298, 301 (2002). 
[44] W. Kwok and Z. Chen, in Proceedings of the Ninth International Meshing Round Table, New 
Orleans LA (2000), pp. 325-333. 



19 



